Antiferromagnetic Zigzag Spin Chain in Magnetic Fields at Finite Temperatures 
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We study thermodynamic behaviors of the antiferromagnetic zigzag spin chain in magnetic fields, 
using the density-matrix renormalization group method for the quantum transfer matrix. We focus 
on the thermodynamics of the system near the critical fields in the ground-state magnetization 
process(M-_H r curve): the saturation field, the lower critical field associated with excitation gap, and 
the field at the middle-field cusp singularity. We calculate magnetization, susceptibility and specific 
heat of the zigzag chain in magnetic fields at finite temperatures, and then discuss how the calculated 
quantities reflect the low-lying excitations of the system related with the critical behaviors in the 
M-H curve. 



PACS numbers: 75.10.Jm, 75.40. Cx, 75.30.Cr 
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Magnetization processes ( M-H curves, where M 
is magnetization and H is magnetic field) of low- 
dimensional antiferromagnetic (AF) quantum spin sys- 
tems have attracted much attention, in accordance with 
remarkable developments in material synthesis tech- 
niques and high-field experiments. Recently, a lot of the- 
oretical and experimental researches have clarified the 
mechanisms of various interesting behaviors of the M-H 
curves at the zero temperature: e.g. the critical phe- 
nomena AAI ~ yJH — H Q associated with the gapped 
excitation (excitation gap cx iJ c )EJ~T] or with the sat* 
urated magnetization tat-,the saturation field H s )p~ u 
magnetization plateaulrfc3 and, the middle-field cusp 
singularity(MFCS)Eirfc-X These are field-induced phase 
transitions of the ground states, reflecting the non-trivial 
energy level structures of the systems. 

Such non-trivial structures of the excitation often 
cause various characteristic behaviors on bulk quantities 
at finite temperatures as well. For instance, thermody- 
namic quantities in magnetic fields ace calculated for lad- 
der systems ,E3 mixed spin systems ,cB etc, where various 
peak-structures of the specific heat are observed. 

In this paper, we study the thermal behaviors of the 
antiferromagnetic zigzag spin chain in magnetic fields, 
which is one of typical systems exhibiting the phase- 
transition behaviors in the M-H curves and actually 
realized as quasi-one-dimensional materials: SrCuC>2 — , 
Cu(ampy)Br2t3 and the organic compound F2PIMNHM. 
However, the finite-temperature properties of the zigzag 
chain in the thermodynamicJimit have not been studied 
quantitatively so muchJ23~E3 since the frustrated inter- 
action makes the ijeliable Quantum Monte Carlo(QMC) 
simulation difficulto. Thus it is a fairly interesting prob- 
lem to investigate the finite-temperature behavior of the 
zigzag chain both from theoretical and experimental view 
points. 

The Hamiltonian of the zigzag chain is given by 



where o, is the S = 1/2 spin operator at i-th site, g is the 
g-factor, and [Ib is the Bohr magneton. We have denoted 
the nearest neighbor coupling as J, the ratio of the next- 
nearest coupling as j(> 0), and the applied field as H. In 
ref. 14, it has been shown that the zigzag chain has fas- 
cinating M-H curves as varying j. Near the saturation 
field, there exists the MFCS for j > 1/4 , and the as- 
sociated two component Tomonaga-Luttinger liquid be- 
havior is observed, where the dispersion curve of the el- 
ementary excitation is the double-well curve. At just 
j = 1/4, the M-H curve behaves as AM ~ {H-JL) 1 ^, 
unlike the usual square-root behavior in j < 1/4.E3 Near 
the zero field, theJi-i? curve in the spin-fluid phase 
j < j fd (= 0.2411)E3 is similar to that of the 5=1/2 
Heisenberg spin chain. For j > jfd, the M-H curve ex- 
hibits the square-root behavior, since the system is in 
the dimerized phase and becomes gapful. Further, for 
j > 0.5, the incommensurate behavior of the ground state 
correlation functior£3~cJ suggests that another MFCS 
appears in the M-H curve near the lower critical field. 
We focus on how the above characteristic behaviors of 
the M-H curves influence the thermodynamic properties 
of the zigzag chain. 

In order to calculate the finite-temperature quantities 
of the zigzag chain, we employ the density matrix renor- 
malization group (DM RiG) metho<£] for the quantum 
transfer matrix(QTM).Eilc3 The remarkable point is that 
the DMRG is free from the negative sign problem and 
thus it is successfully applied to some frustrated spin lad- 
der systems.EilO We calculate the magnetization M, sus- 
ceptibility Xi an d specific heat C at various temperatures 
down to T / J ~ 0.05 in the thermodynamic limit. Then 
we discuss the effects of the non-trivial energy spectrum 
on the obtained quantities. 

This paper is organized as follows. In section 2, we 
describe the DMRG for the QTM in brief. In section 
3, we show the calculated results for the M-H curves of 
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the zigzag chain. In section 4, the thermodynamic prop- 
erties of the zigzag chain near the saturation field are 
discussed in detail. In particular, we consider the effect 
of the MFCS on the specific heat. In section 5, we cal- 
culate the quantities \, C, and C /T near the zero field, 
and discuss the relation to the low-energy excitation of 
the system. Conclusions are summarized in section 6. 

II. NUMERICAL METHODS 

The DMRG method for the QTM is widely used to 
study one-dimensional(lD) quantum spin systems at fi- 
nite temperatures numericallyj2J'E2l since the physical 
quantities can be obtained successfully down to low tem- 
peratures, where the QTM with a sufficient large Trotter 
number should be treated accurately. Further, it should 
be mentioned that the DMRG for the QTM is free from 
the negative sign problem the QMC simulation is suffer- 
ing from. Thus we can say that the DMRG for the QTM 
is one of the most suitable numerical tools to study the 
zigzag chain. 

To see our strategy concretely, we rewrite the Hamil- 
tonian of the zigzag chain with the lattice length 2L into 
a ladder form shown in Fig. 1: 

L 

H = jJ2hi+i, ( 2 ) 

with 

+ j[Sf ■ Sf +1 + §? ■ S? +1 ] + Zeeman terms, (3) 

where A and B are the labels of the lower and upper legs 
respectively. 




FIG. 1. The ladder representation of the zigzag spin chain 

By using the Suzuki- Trotter decomposition, we map 
the zigzag chain at finite temperature into a 2D classical 
system on a checker board lattice.E^I Then the partition 
function of the system can be represented by the QTM 

T (N) T (N)_ 

Z= lim Tr[(T e WT W) L / 2 ], (4) 

N— >oo 

where N is the Trotter number. The QTM is defined by 
the product of the local operator Wi t k originating from 
the local Boltzmann weight in the mapped system: 



t c n) = n w *>*» f ° r *= even ( 5 ) 

2<fc<2JV; feeeven 

T^= 11 W hk , for i = odd (6) 

l<fc<2AT-l; fcGodd 

where i is the index of the spatial direction and k is that 
of the Trotter one. By using the S^-diagonal represen- 
tation |s) of the spin operator S, we have the explicit 
element of the weight: 

/ A.B A,B i w I A,B A,B i _ 
\ s i+l,k s i+l,k+l\ VVi , k \ s i,k s i,k+l) ~ 

I A,B A.B I / f \i A.B A.B v /-n 

(\k+i s i+i,k+i\ e M-£Ki+v\ s i,k s m,fc>' ( 7 ) 

where e = J/(ksTN) and the notation \sfj^) = 
\sf k ,sf k ) is introduced for simplicity. 

In the thermodynamic limit L — » oo, the thermal prop- 
erty of the system is extracted from the largest eigenvalue 
and the corresponding eigenvector of the QTM with the 
sufficiently large Trotter number N . Then we can em- 
ploy the DMRG for the transfer matrix developed by 
T.Nishino.Efl Practically, we regard, for example, \s ik ) 
as a four-state single spin, and perform the DMRG calcu- 
lation with the periodic boundary condition. The detail 
of the algorithm follows Ref.35, where Tc N%> and To are 
renormalized separately. In addition we make the density 
matrix block-idiagonal with help of the conservation law 
for the QTM.EI 

In the DMRG calculation, we often meet undesirable 
complex eigenvalues of the density matrix, because the 
numerical diagonalization yields inaccurate results for an 
asymmetric matrix having degenerate eigenvalues. Thus, 
if the complex number appears in the nearly-degenerate 
eigenvectors of the density matrix, we reorthogonalizc 
the corresponding eigenvectors to be represented by real 
numbers.E3 

We calculate the magnetization M and the internal 
energy E from the obtained eigenvector directly. We fur- 
ther calculate the susceptibility \ an d specific heat C by 
numerical differentiation for M and E respectively. In 
the following, we set ks — I, g^B = 1 and J = 1 for 
simplicity. The calculations were done with e = 0.1 and 
the maximum number of the retained bases m = 88. We 
have confirmed that the computed data converged with 
respect to m and e. 



III. NUMERICAL RESULTS FOR THE M-H 
CURVES 



In Fig. |^, we show the calculated M-H curves of the 
zigzag chain at finite temperatures for j = 0.2, 0.5, and 
0.6. We also show the zero-temperature M-H curves 
calculated by the product-wavefunction renormalization 
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group method (PWFRG)il for comparison. At zero tem- 
perature, the M-H curve for j = 0.2 has no anomaly in 
the middle field region. On the other hand, the M-H 
curve for j = 0.5 has one MFCS at H, 



1.9. We fur- 
ther find that the M-H curve for j = 0.6 has two MFCS 
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FIG. 2. the M-H curves of the zigzag chains at finite tem- 
peratures, (a) j = 0.2, (b) j = 0.5, and (c) j — 0.6. Inset: 
the differentiation of the M-H curve (dM/dH curve). 

At finite temperature, the singularities of the ground 
state M-H curve are generally rounded by the thermal 
excitation. Indeed the M-H curves show no singular- 
ity in the high temperature region. As temperature de- 
creases, however, the quantum effects appear. We can 
see that the temperature dependence of the M-H curve 
for j = 0.5 is enhanced near the cusp field H cusp ~ 1.9, 



comparing with that for j — 0.2. The differentiation 
curve(dM/dH curve) of the M-H curve shows the MFCS 
effect more clearly (insets of Fig. ||). In the dM/dH curve 
for j — 0.5 at T — 0.05, a shoulder can be observed at 
H = 2.0. For j = 0.6 we can see the winding structure 
in the M-H curve at T = 0.05, corresponding to the 
higher MFCS and the lower MFCS. The dM/dH curve 
also displays the double cusp structure of M-H curve at 
low temperature; the shoulder at the higher MFCS field 
becomes more clearly and the peak near the lower MFCS 
field does sharp. 



IV. NEAR THE SATURATION FIELD 



In this section, we consider the thermodynamic behav- 
iors of the zigzag chain near the saturation field, where 
we can take a down-spin-particle picture of the elemen- 
tary excitation.Ej The saturation field H s is given by 2 
for j < 1/4 and 1 + 2j + l/(8j) for j > 1/4. Near H s , a 
down spin in the saturated (all up) state can be regarded 
as a spinless fermion, which is the effective low-energy 
limit of the (5-function Bose gas modeD'H, and then the 
energy dispersion of the particle is calculated to be 



u(k) = cos A: — 1 + j(cos2fc — 1). 



(8) 



This one-particle dispersion curve fully characterizes the 
qualitative property of the M-H curve near H s . For 
j < 1/4, u(k) has a single minimum at k — n, while, 
for j > 1/4, Lu(k) has a local maximum|-at-|fc = n and 
two minima at k = n ± cos _1 (l/(4j)).E3'til Thus, for 
< j < 1/4, the M-H curve is smooth in the whole 
field range < H < H s . While, for j > 1/4, " van Hove 
singularity" corresponding to the local maximum of u(k) 
gives jttje MFCS in the higher-field region of the M-H 
curve! 
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FIG. 3. The specific heat C near H = H s . (a) j = 0.2 and 
(b) j = 0.5. Inset of figure (a): C/T 1/4 plot for j = 0.1, 0.2 
and j = 1/4 at H = H 3 . Inset of figure (b): magnification in 
low-temperature region at the slightly upper fields than H a . 

In Fig. ||, we show the specific heat C at various 
applied fields H above H s . In the case of j = 0.2, 
clear peaks of C can be seen in the low-temperature 
region(T < 0.3), except for H = H s . Then we can find 
that the peak position Tpoak, which is defined by the 
temperature at the top of the low-temperature peak, is 
approximately proportional to the energy gap \H — H s \. 
Hence, the low-temperature peaks represent the diver- 
gence of the density of state at the bottom of the disper- 
sion curve (||). 

At just H — H s , the peak disappears and, then, the 
temperature dependence of the specific heat C{H = H s ) 
becomes the power-law behavior: T 1 /" in T <C l(see the 
inset of Fig. |(a) ). This power a is determined by the 
shape of the dispersion curve in k — > ir limit, where we 
have ui(k) ~ (k — 7r) Q . For j < 1/4, in strict sense, a = 2 
should be obtained in T — > limit, since Eq. (H) has a 
quadratic term. However, the coefficient of the quadratic 
term decreases and the quartic term becomes dominant, 
as j increases to j — 1/4. Thus we see the cross-over 
behavior from a = 4 to 2 as shown in the inset of Fig. 
3(a). On the other hand, C/T 1 / 4 plot for j = 1/4 is con- 
sistent with a = 4, where the quadratic term vanishes 
completely and Eq. (^) becomes uj(k) ~ (k — tt) 4 . The 
above power-law behavior of C is essentially the same as 
the Fermi liquid one, which is verified for the Heisenberg 
chain(j = 0) with a = 2Ej. 

For j — 0.5, the double-well structure of Eq. (||) in- 
duces interesting properties on the specific heat C near 
H s . Indeed, the low-temperature peak exists even at 
H = H s (= 2.25), and a weak shoulder structure can 
be observed at the slightly upper fields than H s (see the 
Inset of Fig. ||(b) ). In addition, we can see that the 
peak-temperature T poa k is approximately proportional to 
H — H cusp when \H S — H\ <C 1. Therefore the strong peak 
represents the "van Hove singularity" corresponding to 
the MFCS, and the weak shoulder structure comes from 
the bottom of the dispersion curve. 




FIG. 4. The specific heat C near H — H cusp . For com- 
parison, we show that for j = 0.2 at H — 1.6. 

In order to illustrate the effects of the MFCS more 
clearly, we have calculated the specific heat C near the 
cusp field -f/cusp — 1-9. In Fig. |], we show the calculated 
result for j = 0.5 at H = 1.85(= H s - 0.4), which is 
slightly below H cusp . We also show that for j = 0.2 at 
H = 1.6(= H s — 0.4) for comparison. For j = 0.5, we 
can find that the specific heat increases sharply at very 
low temperature, and has a shoulder at T ~ 0.1. On the 
other hand, for j = 0.2, we find a peak near T ~ 0.1 only. 
Thus, the sharp increase of the specific heat at very low 
temperature for j = 0.5 comes from the "van Hove sin- 
gularity" of uj(k), and the shoulder is attributed to the 
bottom of Lu(k). In addition, we note that the bottom of 
u)(k) gives weak contribution to the specific heat. 

From the above results, we can conclude that the band 
edge singularities of ui(k) at k = n or k = 7r±cos -1 (l/4j) 
explain correctly the specific heat behaviors both near H s 
and i?cus P - 



V. NEAR THE ZERO APPLIED FIELD 



At the zero magnetic field, varicjas [experiments have 
been done for the actual materialstlrtj, where the ob- 
served susceptibility is shifted from that of the pure 
5 = 1/2 Heisenberg model. This shift of the suscepti- 
bility is attributed to the deformation of the low-lying 
excitation, which is induced by the frustration effect due 
to the next-nearest coupling. In this section, we calcu- 
late the thermodynamic quantities near the zero field, 
which play a crucial role to determine the exchange cou- 
pling constant of the materials, and discuss the relation 
between the obtained quantities and the low-lying exci- 
tation structures. 

In Fig. ||, we show the specific heat C, the quan- 
tity C/T, and the susceptibility \ a t the zero field for 
0.1 < j < 0.6. 
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FIG. 5. Thermodynamic quantities at zero magnetic field 
(a) The specific heat C, (b) C/T, and (c) The susceptibility 
X- 

In j < jid(= 0.2411), the zigzag chain is in the gapless 
spin-fluid phase at the zero field and the M-H curves 
do not show any singularity as in Fig. ||(a). Thus the 
calculated quantities for j — 0.1 and j — 0.2 are qual- 
itatively the same as that of the 5=1/2 Hcisenberg 
chain(j = 0), though the broad peaks of x and C shift 
to the low-temperature side. As shown in Fig. ||(a) 
and (b), x and C/T approaches to the constant value 
in the zero-temperature limit. We extrapolate the coef- 
ficient of the linear specific heat 7 = linvr-^o C /T by a 
polynomial fit of C/T = 7 + a{T + C12T 2 in the range 
0.06 < T < 0.12, and obtain 7 = 0.75 for j = 0.1 
and 7 = 0.85 for j — 0.2. We also compute the zero- 



temperature susceptibility xo by numerical differentia- 
tion of the PWFRG results with AH = 0.005: we obtain 
Xo = 0.118 for j = 0.1 and xo = 0.131 for j = 0.2. Here 
we note that the same manner of the extrapolation at the 
Heisenberg point (j = 0) yields 7 = 0.66 and xo — 0.108, 
which agree with the Bethe ansatz values. According to 
the conformal field theory (CFT), the coefficient of the 
linear specific heat and the ground state susceptibility 
Xo are given by 7 = nc/(3v) and xo = 1/(2ttv) respec- 
tively, where c is the central charge and v is the spin 
wave velocity. The calculated values of 7 and xo are al- 
most consistent with the CFT pxediciipn with c = 1 , but 
the strong log-correction for Xo"'i— MSI prevents us from 
verifying the CFT relation precisely. 



For j > jfd the zigzag spin chain is gapful. How- 
ever, the [Magnitude of the gap is quite small in jfd < 
j < 0.4ESE3 Thus the specific heat C and the suscepti- 
bility x f° r j = 0.3 are similar to those in j < j^. For 
j = 0.4 and 0.5, x and C catch the effect of the gap: 
the exponential decay of the susceptibility can be seen 
in the low temperature. For j = 0.6, the energy gap 
becomes large, so that x an( i C show clear exponential 
decay in low temperature region. In the specific heat for 
j = 0.6, the peak at T ~ 0.15 seems to be distinguish- 
able from the broad peak around T ~ 0.5 (see Fig. ||(a)). 
As shown in Fig. |^(c), the M-H curve of j — 0.6 has 
the lower MFCS, and the dispersion curve of the excita- 
tion is expected to be the double-well structure, which 
is supported by the incommensurability of the ground- 
state correlation function.E3~E3 The double- well structure 
of the dispersion curve possibly explains the remarkable 
shape change of the specific heat for j = 0.6. 



Further, we consider the zigzag chain in a small mag- 
netic field, where the excitation related with the lower 
MFCS can be observed more directly. In Fig. ^, we show 
the specific heat C and the susceptibility x f° r 3 = 0.2, 
0.5,and 0.6 at the fixed field H = 0.7. For j = 0.2, the 
properties of C and x is essentially the same as those at 
the zero field. For j = 0.5, on the other hand, C and x 
are enhanced in the low temperature region, because of 
the excitation gap. For j = 0.6, where the lower MFCS 
appearing in the M-H curve, we can see the susceptibil- 
ity has very sharp peak at very low temperature. The 
specific heat also grows sharply at very low temperature 
as well, which is quite similar to that of the higher field 
MFCS. These outstanding peaks in x and C support that 
the dispersion curve of the low-lying excitation has such a 
structure as the double-well one. However, the evidence 
for the double- well structure can not be detected directly 
within the present calculation, since the field at the lower 
MFCS and the lower critical field of the excitation gap 
are very close as shown in Fig. ||(c) . 
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FIG. 6. The specific heat C and the susceptibility \ m the 
magnetic field H = 0.7. 

Here, we make a comment on the relevance of present 
calculations with the experimental results. For exam- 
ple, the susceptibility of the compound Cu(ampy)Br2 was 
measured by Kikuchi, and then the_ratio of the exchange 
coupling is estimated as j = 0.2.E3 The M-H curve of 
this compound is also observed to have no anomaly in the 
middle field region, which is consistent with the present 
calculation. Moreover, Hosokoshi, et al have found that 
the susceptibility of the organic compound F2PIMNH 
shows .-the clear exponential decay in low temperature 
regionta, which imply j > 0.4 compared with the calcu- 
lated susceptibility in Fig. ||-(c). Then it is suggested 
that the M-H curve of F 2 PIMNH has the MFCS. How- 
ever the high field experiment has not yet been done for 
this compound. 



VI. SUMMARY 



In this paper we have quantitatively studied the ther- 
modynamic behaviors of the antiferromagnetic zigzag 
spin chain in magnetic fields. We have calculated the 
magnetization processes {M-H curves) at finite tem- 
peratures by using the density-matrix renormalization 
group (DMRG) method for the quantum transfer ma- 
trix. The zero-temperature M-H curves were also cal- 
culated by the product wave-function renormalization 
group method. We have shown that the zero-temperature 
M-H curve has zero, one, and two middle field cusp sin- 
gularity(MFCS) for j = 0.2, 0.5 and 0.6, respectively. 
The thermal effect on the M-H curve was discussed as 
well. 

We have further investigated the bulk physical quan- 
tities at finite temperatures in terms of the non-trivial 
energy-level structures responsible for the singularities 
in the zero-temperature M-H curve. Near the satura- 
tion field, we have considered the band edge singularities 
m i < 1/4, where the Fermi liquid behavior is observed. 
For j > 1/4, we have shown that the double-minimum 
shape of the dispersion curve correctly explains the peak 
structures of the specific heat. 



Near the zero field, we have considered the suscepti- 
bility x an< 4 specific heat C. In the gapless region, i.e. 
j < 0.2411, we have shown x an( i C exhibiting typical 
gapless behaviors and checked the CFT relation between 
the zero-temperature susceptibility and the coefficient of 
linear specific heat. For j > 0.2411, the exponential de- 
cay of x an d C has been observed in low temperatures. 
Moreover, we have seen that another peak of C is in- 
duced for j = 0.6, where the double- well structure of 
the dispersion curve is expected to accompany with the 
lower MFCS. We have also calculated the susceptibility 
X and specific heat C in the small field H — 0.7. The 
highly enhanced peaks for \ and C of j = 0.6 support 
the double-well structure of the dispersion curve. To un- 
derstand the microscopic view about such low-lying exci- 
tation connected to the conformation mechanism of the 
lower MFCS is a remaining problem. 

In the connection to the experiments, the M-H curve 
of F2PIMNH is interesting, since the appearance of the 
MFCS can be predicted by the present calculation. In acL 
dition, the crystal structure of another zigzag materials^ 
suggests that they may have large next-nearest coupling 
j, by taking into account the orbital symmetry of the 
atoms concerned with the super-exchange interaction. In 
analyzing experimental data for the above zigzag mate- 
rials, we believe that the present results is of use. 
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